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ABSTRACT 

An approximate temperature of cr^'stallization of rocks and im- 
portant information on the thermodynamic nature of a crystalline solu- 
tion, coexisting with orthopyroxene, can be obtained by using the experi- 
mental Mg^^-Fe^* intra-crystalline distribution data in orthopyroxenes 
at various temperatures and the Mg^*-Fe^*^ inter- crystalline distribu- 
tion data in natural assemblage. Ca-pyroxene Ca (Mg, Fe) Si^O^ with 
the M2 site nearly filled by Ca, is a 'simple mixture' (Guggenheim, 
1967) of Mg^* and Fe^^ components with an energy constant w approxi- 
mately 1100 cal/mole at 600°C. The estimated temperature of ion- 
exchange equilibrium between coexisting cummingtonite and orthopyro- 
xene in Biwabik Iron Formation, Minnesota is ~610°C. Cummingtonite 
also approximates to a 'simple mixture' with w as 700 cal/mole at 


600°C. 


INTRODUCTION 


The distribution of and between two cc existing ferromagnesian 

minerals such as ortho- and calcic-pyroxene by Kretz (1963) has been studied 
on the assumption that the Fe and Mg components mix Ideally in each of the 
minerals. Mueller (1961, 1962) considered the possibility of intra-crystalline 
partitioning of the cations among the different structural sites in a mineral and 
for cummingtonite found it desirable to assume a *r"';ular solution’ model. 

Matsui and Banno (1965), Banno and Matsui (1966) and Grover and Orville (1969) 
assumed ideal mixing of Fe*^ and Mg 2+ on each site in orthopyroxene and 
interpreted the inter-crystalline Mg ^^-Fe distribution between coexisting 
orthc- and ca-pyroxene on that basis. 

Saxena and Ghose (1970a) have shown that the mixing of Fe^^ and Mg^^ on 
Ml and M2 sites is better approximated by a ’simple mixture’ model (Guggenheim, 
1967) for each site individually. The orthopyroxene crystalline solution as a 
whole is found to be non- ideal asymmetric particularly at 500 and 600° C (Saxena 
and Ghose, 1970b). By using the intra-crystalline Mg^^-Fe^* distribution data 
in orihopyroxene and the information on its solution behavior, the inter-crystalline 
Mg^^-Fe^^ distribution between orthopyroxene and another coexisting mineral 
CEin be interpreted more reliably. The study of Mg^*-Fe^* distribution between 
orthopyroxene and another mineral yields Important information on the ternpe *a- 
ture of crystallization of the host rock and on the thermodynamic nature of mixing 
of the Mg and Fe components in the other mineral. 
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The following notations and symbols are used; 

Opx - orthopyroxene, Cpx- calcic pyroxene 
Cum - cummlngtonlte 

- equilibrium constant for a reaction (a) 

- mole fraction Fe^*/(Fe2^ + Mg^^) in a phsise a as on a site a in a phase 
a® - actlvltj^ of a component A in a phase a or on a site a ; in the latter case 

they are referred to as 'partial* activity 
w - an energy constant for 'simple mixture' (see Guggenheim, 1967) 

AG° - standard free energy change for a reaction (a). 

COEXISTING ORTHO- AND CALCIC -PYROXENE 
The Ion-exchange Reactions 

The ion-exchange reactions between orthopyroxene (Fe, Mg) Si O with Ml 

2 2 6 

and M2 sites and Ca- pyroxene C a (Mg, Ft) Si O with M2 site completely occupied 

2 6 

by Ca are: 


Mg (Ml) + Fe (Ml) 

^Fe (Ml) + Mg (Ml) 

(1) 

Opx 

Cpx 

Opx 

Cpx 


Mg (M2) + Fe (Ml) 

r^Fe (M2) + Mg (Ml) 

(2) 

Opx 

Cpx 

Opx 

Cpx 


Mg (Ml) + Fe (M2) 

(Ml) +Mg (M2) 

(3) 

Opx 

Opx 

Opx 

Opx 
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The equilibrium constants for these reactions are: 




gMl gCpx 

Fe 

gMl flCpx 

Mk Fo 


*2 = 




*3 = 


-Ml -M2 
^Fe ^Mg 

gMt ^M2 
Mr Fe 


where represents the partial activity of referred to the Ml site In 
orthopyroxeiie and similarly ajJJ , ajlj ai^d a[[“ represent other partial activities, 
a^* or a^P* are activities of the Fe and Mg components In Ca-pyroxene. In 
the present case, they are the same as the partial activity of Fe or Mg on the 
Ml site In Ca-pyroxene. The Inter-crystalline Ion-exchange reaction between 
(Mg, Fe) SiO, and Ca(Mg, Fe) SLO. Is 

3 2 0 

Mg (Ml) + 1/2 Fe (Ml + M2) ^ Fe (Ml) + 1/2 Mg (Ml + M2) W 

Cpx Opx Cpx Opx 

The equilibrium constant for (4) Is 


_ a°;» 


a^x (aMi 
Fe ' Mg 



1/2 


aSr (a 


Ml 

Fe 


^«2)K2 


( 5 ) 


Transforming K In terms of K, and K- i we have 

4 1 ' 




and are Interrelated as 


*3=*;. 




O) 
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As discussed by Kretz (1963) and Grover and Orville (1969), the effect of pressure 
on these equilibrium constants Is small. Therefore, we shall be mainly concerned 
with the temperature effect. 

Distribution of and Fe between Ca-pyroxene and Ml Site In 

Coexisting Orthopyroxene 

The Intra-crustalllne Mg^^-Fe^* distribution data between Ml and M2 sites 
In orthopyroxene (Saxena and Ghose, 1970a) can be used to obtain the Inter- 
crystalllne distribution between Ml or M2 site In orthopyroxene and Ml site In 
rllnopyroxene If the equilibrium temperature of the Inter-crystalline Ion-exchange 
Is known. In absence of any experimental data on the Mg^"*^-Fe^* distribution 
between ortho- and C a- pyroxenes In the temperature range of 500-800®C, the 
data available from natural assemblages are used. This Involves both uncertainty 
In the temperature of orystalllr.atlon as well as the non-stolchiometry of Ca In 
the M2 site In cllnopyroxene. Both the pyroxenes In natural assemblages also 
contain small amounts of Al^* , Fe^"*^, Tl^* and Mn^* . All these factors shall 
affect the quantitative results obtained by using the chemical data I'rom natural 
assemblages. 

The data on the distribution of Fe^^and Mg^^ between coexisting ortho- and 
Ca-pyroxene in six rocks of granutite facies, chosen arbitrarily from the data 
In literature, are presented In Table 1. Assuming that the temperature of Inter- 
crystalllne equilibrium was 600®C In these rocks, the data on In orthopyroxene 
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correspoLcU.'ig to each composition are obtained from Fig. 1 (reproduced here 
from Saxena and Ghose, 1970a, Fig. 1). Similarly, the data on vJJ in ortho- 
pyroxene at 700‘’C are obtained and listed in Table 1. The distribution coefficient 
K\ is; 



1 - yCp* 

^ ''Fe 


yCpx 

“'re 


VMI 

Fe 




( 8 ) 


would correspond to the equilibrium constant if both Ca-pyroxene (Ml 
site) and the Ml site in orthopyroxene were ideal. The values of K[ in Table 1 
show that one or both of the phases may be non-ideal. As the Ca-pyroxene in 
the six samples do not contain Ca^*; (Mg^“^ + Fe^*) as 1:1, the values may 
change as a function of Ca concentra*’.on. A plot of concentration of Ca^*^ in Ca- 
pyroxene against does not show any significant correlation between the two. 
Fig. 2, however, shows that there is a definite correlation between K'^ and . 


Although Ca^* does not completely fill M2 site In the slxCa-pyroxenes, a 
one site quasi-binary model for Ca-pyroxene may still be assumed. The absence 
of any correlation between XJ and concentration of Ca^^ in the pyroxene favors 
this assumption. 


Since mixing on the Ml site in orthopyroxene is non-ideal (Saxena and 
Ghose, 1970a), the K[ values can be improved by considering another equilibrium 
constant K[‘ as 
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(9) 


K 


n 

1 “ 


> - xfV 


e 


Ml 

Pe 



These values, also listed in Table 1, continue to be a function of , Therefore, 
the mixing of Fe- and Mg-compponents in Ca-pyroxene should be non-ideal to a 
certain extent. 


Estimation of Temperature of Crystallization 

Figure 3 shows a plot of in orthopyroxene against XJp* . The opon 
circles are based on the assumption that the temperature of inter-crystalline 
ion-exchange equilibrium was 600° C. Similarly, the solid circles represent 
700°C. The Ml site in orthopyroxene is a 'simple mixture' (Saxena and Ghose, 
1970a) of Fe^* and Mg^^. Assuming that the one site quasi-binary Ca-pyroxene 
is also a 'simple mixture', the log of the equilibrium constant for (1) is 


11/Ml u/Cpx 

in ^ = In a:; + _ (1 . 2 XMJ ) - ( 1 - 2 ) 


RT 


( 10 ) 


Equation (10), in an exponential form, was derived by Mueller (1964) and later 
used by Saxena (1969) in the logarithmic form. In Fig. 3, the distribution curves 
are drawn using a least squares fit to the distribution data in Tabie 1 according 
to equation (10). Table 2 shows the resulting constants. From Saxena an*' Ghose 
(1970a), the experimental values of w**^/RT are 1.294 an 007 at 600 and 
700°C respectively. The difference between the experimental /RT and the 
calculated w**V/?r in Table 2 is due to the actual temperature of crystallization 
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being different from either 600 or 700®C. By plotting the experimental w**^/RT and 
the calculated /RT In Table 2 against temperature, the temperature of 
crystallization can be determined graphically from the Intersection of the straight 
lines, one Joining the calculated /RT and the other joining the experimental 
w^^/RT. This method of temperature estimation is quite geneial and c^n be 
applied to a group of coexisting pyroxenes with sufficient degree of Mg;Fe varia- 
tion In the range of temperatures for which the Intra-crystalllne values are 
available. The present data yield 620®C as the temperature of Inter-crystalllne 
Ion-exchange equilibrium (4). Corresponding to this temperature Is 

0.625 and Is 0.85. 

Mg Partitioning In Metamorphlc and Igneous Pyroxenes: 

Assuming that both orthopyroxene and Ca-oyroxene ap’^roach ideal solution 
closely at high temperatures corresponding to igneous rocks (1150®' >, can 
be taken as 0.73 as obtained by Kretz (1963). w^^^/RT at 600°C Is then 0.650 
(or w is 1100 cal/mole) and the activity-composition relation (Fig. 4) In Ca- 
pyroxene Ca (Mg, Fe) Sl^Og Is: 

0-650 (1 - 

where afp* refers to the activity of Fe?* in the pyroxene crystalline solution 
and should be similar In magnitude to the activity of the end-member CaFejSi^Og. 
The equation (11) Is for the ’simple mixture' as defined by Guggenheim (1967). 
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By taking a» 0.85 at 600®C md combining it with as 0.2516 (from the 
experimental results in >:;’axena and Ghose, 1970a), is 3.38. is calculated 
to be 0.59 and AO® is 915 cal/mole. The distribution curve at 600® C tor the 
inter-crystalline Mg ^*-Fe distribution according to these data is shown in 
Figure 5. The ideal distribution curve with as 0.73 is also plotted. ^ 

It was noted by Saxena and Ghose (1970a) that the intracrystalline ion- 
exchange free energy AC® does not change significantly wiia changing tempera- 
ture. It is probable that AG® also does not vary much with temperature. If 
we assume 1150®C for as 0.73, AG® is 890 cal/mole v'hich is of the same 
order of magnitude as AG® at 600®C The different distribution isothernu 
in the range of 600-1200 ’C shall be quite closely spaced. As a result, they 
cannot yield the temperarure estimates very accurately. 

COEXISTING (>PTHO PYROXENE AND CUMMINGTONITE 
Estimation of Temperature 

The data on the distribution of Mg^"^ and Fe^* between coexisting ortho- 
pyroxene and cummingtonite in twelve metamorphic samples from Biwabik 
Iron Formation, Minnesota (Bonnichsen, 1969) are listed in Table 2. The ion- 
exchange r eaction between cummingtonite and Ml site in orthopyrovvene is written 
as 

^Fe-(all sites) + Mg-Ml^iMg - (all sites) + Fe - Ml 
Cum Opx Cum Opx 


8 



Disregarding the different sites in cummingtonite and assuming that the crystal 
as a whole is a 'simple mixture', the equilibrium constant for (12) is 

InKjj =lnjCp + _ (1-2X“J)-— (l-2XCp (13) 

Assuming that the temperature of equilibrium in (12) is 500, 600 and 700°C, a 
least squares fit according to (6) for the data in Table 2 at each temperature 
is obtained. The results are shown in Table 3. The values of calculated /rt 
in Table 3 are plotted against temperature and so also the experimental values 
of in orthopyroxene from Saxena and Ghose (1970a). A crossing of the 

two curves gives the temperature of inter-crystalline equilibrium (12) as 610°C. 

The Mg^^-Fe^^ Distribution between Coexisting Orthopyroxene and 
Cummingtonite at 600° C 

At 600 °C the approximate value of wCum is 0.42. The activity of Fe^"^ in 
cummingtonite at 600 °C is: 

^-42 (1 - (14) 

The value of the energy constant w is 728 cal/mole at 600°C. The activity of 
Fe component in orthopyroxene may be calculated by using the anal 3 rtical expres- 
sion for the activity-coefficient (Saxena and Ghose, 1970b) 

+ 0-118 (3 X,. - + 0.190 (X,. - X„,) (5 
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or {rtphlcally from Fig. 3 la Saxena and Ghose (1970b). The ioa-exchange 
raactlon batweeo coexlaUng curomlngtonite and orthopyroxene is: 


- Mg - Cua ♦ Fe - Opx 
7 


Mg - Opx 4 1 Fe • Cum 


( 16 ) 


The equi librium constant for (16) is 



Cwa jkOpo 

r. 


aCu. 

r» F» 


(17) 


TU>lo 4 shows Uu value of K in six samples, selected to represent a wide 

compoeition range. The values of using as are also listed 

for comps, son. The values of g are more uniform but they still vary with 

1 0 

or It is probable that cummingtonite, similar to orthopyroxene 

(Saxena and Ghose, 1970b) is asymmetric crystalline solution and not the 'simple 
mixture' assumed in the preceding section. However, the deviations from a 
'simple mixture' model is probably not great. The at 600^C, assuming 
as 0.79, is 400 cal/mole. 

DISCISSION 

The availability of the intra-crystalline Mg^^ -Fe^*^ distribution data in 
orthopyroxene at various temperatures when used with the Mg^'^-Fe^'^ distribu- 
tion data for the coexisting ortho- and Ca-pyroxene in natural assemblages 
makes it possible to interpret such inter-crystalline Mg^^-Fe^^ distributions 
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without making various assumptions which were necessarily used by previous 
workers. In particular, it is not necessary to assume 1) the ideal mixing of 
Mg^* and Fe^^ on individual sites in a p 5 rroxeae and 2) the equivalence of 
Ml site in Ca-pyroxene to the Ml site in orthopyroxene. The second assumption 
was introduced by Mueller (1962) and was used along with the first assumption 
by Banno and Matsui (1965) and Grover and Orville (1969). Although these as- 
sumptions are no longer required in the present work, a binaiy solution model 
for a multi-component phase has still to be used. Ca-pyroxene should be treated 
as a phase with three major components (Fe^"^ , Mg^* and ) distributed 
between two non-equivalent sites. The solution model for cummingtonite, a 
quasi-binary four site phase, is not yet definite. 
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TABLE 1 


Distribution of and Mg between Coexisting Pyroxenes 


Sample 



In Orthopyroxene 


-Ml 



H-4 

0.650 

0.773 

0.645 

0.759 

0.608 

0.978 

0.672 

115 

0.507 

0.632 

0.405 

0.640 

0.735 

0.661 

0.846 

4642 A 

0.429 

0.571 

0.310 

0.574 

0.781 

0.598 

0.978 

Arendalite 

0.365 

0.500 

0.219 

0.482 

0.83:t 

0.488 

1.009 

2941 

0.321 

0.457 

0.196 

0.425 

0.855 

0.483 

1.031 

2270 

0.249 

1 

0.384 

0.129 

1 

0,317 

0.900 

0.446 

1.062 

700®C 

H4 

0.650 

0.773 

0.660 

0.741 

0.527 

1.045 

0.757 

115 

0.507 

0.632 

0.435 

0.599 

0.683 

0.748 

0.852 

4642A 

0.429 

0.571 

0.350 

0.535 

0.735 

0.717 

0.964 

Arendalite 

0.365 

0.500 

0.260 

0.450 

0.782 

0.610 

1.001 

2941 

0.321 

0.457 

0.225 

0.412 

0.815 

0.610 

1.068 

2270 

L 

0.249 

0.384 

0.156 

0.319 

0.865 

0.557 

1.114 


Samples: H^-Metamorphic iron formation; northern Quebec; Kronclc (1961). 

115,4642 A, 2941 and 2270-Charnockites from India; Howie (1955). 

Arendalita-From the type locality Arendal, southern Norway; collected by Saxena. 

X = Fe^V(Fe^^ +Mg^^, 3pJ partial activity referring to activity of Fe^^on Ml 
orthopyroxei.e; the values are obtained from Fig. 1 of Saxena and Chose (1970a). 




-Ml 

r"- ff* 

i-xc;» 

, - X«J 


^Mg 
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TABLE 2 


Fe^^-Mg^^ Distribution Data in Coexisting Cummlngtonlte and Orthopyroxene 


Sample Nos. 

yCum 

^Fe 

vOpx 

in Opx 

500°C 

600°C 

700®C 

026 

0.666 

0.736 

0.540 

0.565 

0.580 • 

073 

0.573 

0.631 

0.365 

0.405 

0.435 

082 

0.722 

0.771 

0.635 

0.650 

0.660 

087 

0.732 

0.754 

0.585 

0.600 

0.615 

088 

0.609 

0.637 

0.370 

0.410 

0.440 

122 

0.747 

0.787 

0.650 

0.660 

0.670 

131 

0.274 

0.310 

0.045 

0.085 

0.115 

175 

0.676 

0.750 

0.585 

0.600 

0.615 

234 

0.743 

0.791 

0.670 

0.680 

0.690 

241 

0.659 

0.722 

0.535 

0.555 

0.570 

266 

0.594 

0.634 

0.365 

0.405 

0.435 

304 

0.353 

0.416 

..085 

0.1^0 

0.185 


Data on coexisting minerals from Bonnichsen (1969). Data on XpJ in orthropyxene from 
Fig. 1 (Saxena and Chose, 1970a). 
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TABLE 3 


Equilibrium Constant and w for the Data in Table 2 


\ 


Using a Least Squares Program on Equation 



Xj2 

^Cum /p’f 

w^^/RT (opx) 

500 

0.36 

-1.00 

0.58 

600 

0.61 

0.42 

1.19 

700 

0.78 

1.22 

1.63 


TABLE 4 


Calculation of Equilibrium Constants 


S. No. 

yCum 

Fe 

-Cum 

^Fe 

-Cum 

3^. 


It. 

B 

131 

0.274 

0.310 

.341 

.749 

.420 

.725 

.78 

.65 

304 

0.353 

0.416 

.421 

.681 

.530 

.640 

.75 

.66 

073 

0.573 

0.631 

.620 

.490 

.740 

.430 

.73 

.78 

026 

0.666 

0.736 

.698 

.402 

.810 

.360 

.77 

.88 

082 

0.722 

0.771 

.746 

.346 

.825 

.330 

.86 

1.04 

234 

0.743 

0.787 

.764 

.324 

.840 

.310 

.87 

1.07 


★ 


-Cum 
^ F. 

-Cum 

^8 


gOpx 

Mb y , 


F* 


F* Mg 


*Sr 
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FIGURE CAPTIONS 


Figure 1. Distribution of Mg and Fe^"^ between Ml site and M*’ site in ortho- 
pyroxene (from Saxena and Ghose, 1970a). The mole fractions and 
X^l are determined by Mossbauer resonance spectroscopy on heated 
orthopyroxenes . 


l(\ - X**^ \ 

Figure 2. The distribution coefficient Kl \ —• — — — ) plotted against Xp^* 

/\ x^p* 1-^;/ 

to show the non-ideal relationship. 

Figure 3. The distribution of Fe^^ and Mg between Ml in orthopyroxene and 
clinopyroxene. The curv'es are least squares fit according to the 
equation employing simple mixture model of Guggenheim (see Saxena 
and Ghose, 1970). The dashed curve uses the data on X“^ in ortho- 
pyroxene at 600°C and the solid curve at 700°C. 


Figure 4. Inferred activity-composition relation in calcic pyroxene CafMg, Fe) 
SijOg at 600°C. 

Figure 5. Distribution of Fe^^ and Mg between ortho- and calcic-pyroxene. The 
aata are from Kretz (la63). Open circles - metamorphic, crosses — 
igneous. The 600°C isotherm is based on activity-composition rela- 
tions in clinopyroxene as shown in Fig. 3 and in orthopyroxene as 
listed in Table 1 (Saxena and Ghose, 1970b) is 0.59. The 1150°C 
isotherm is based on ideal solution model with K, as 0.73. 
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Fig;ure 1. Distribution of Mg ^ and Fe‘ between Ml site and M2 site in ortho 
pyroxene (from Saxena and Chose, 1970a). The mole fractions Xj,” i 
X^l are determined by Messbauer resonance spectroscopy on hentod 
orthopyroxenes. 








0.2 0.6 1.0 


V Ml . _ 

Xpg inOpx 

Figure 3, The distribution of Fe^'*’ and Mg between Ml in orthopyroxene and 
cli nopyroxene. The curves are least squares fit according to the 
equation employing simple mixture model of Guggenheim (see Saxena 
and Ghose, 1970). The dashed curve uses the data on in ortho- 
pyroxene at 600®C and the solid curve at 700°C. 








